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We present an experimental investigation of ultrafast-laser driven Rabi oscillations of atomic 
rubidium. Since the broadband spectrum of an ultrafast laser pulse simultaneously couples all 
the electronic hyperfine transitions between the excited and ground states, the complex excitation 
linkages involved with the D1 or D2 transition are energy degenerate. Here, by applying the Morris- 
Shore transformation, it is shown that this multi-state system is reduced to a set of independent 
two-state systems and dark states. In experiments performed by ultrafast laser interactions of 
atomic rubidium in the strong interaction regime, we demonstrate that the ultrafast dynamics of 
the considered multi-state system is governed by a sum of at most two decoupled Rabi oscillations 
when this system interacts with ultrafast laser pulses of any polarization state. We further show 
the implication of this result to possible controls of photo-electron polarizations. 

PACS numbers: 32.80.Qk, 78.47.jh, 42.65.Re 


I. INTRODUCTION 


The Morris-Shore (MS) transformation [T] offers a the¬ 
oretical framework to minimize the complexity in a va¬ 
riety of coherent excitation coupling linkages of a N- 
state quantum system {2J. Often used in modeling 
laser-induced excitation of atoms, the MS transforma¬ 
tion replaces the complex system by a set of indepen¬ 
dent two-state systems, decoupled from each other, and 
unlinked states. For example, when a degenerate two- 
state atom possessing angular momenta J and J ', respec¬ 
tively, for the ground and excited levels, interacts with 
a laser field, the excitation degeneracy in the magnetic 
sub-levels results in various coupling linkages amongst 
N = 2(J + J' + 1) levels, characterized by polarization 
states of the laser field. However, a consequence of the 
MS-transformation decoupling is that the N x N Hamil¬ 
tonian that describes the coherent excitation is reduced 
to a set of 2 x 2 sub-Hamiltonians plus a diagonal matrix 
representing uncoupled dark or spectator states. So, the 
dynamics of such system can be described simply by a su¬ 
perposition of independent two-state dynamics, each of 
which undergoes a Rabi oscillation between a pair of new 
ground and excited states defined in a factorized Hilbert 
sub-space. 

The procedure of the MS transformation can be sum¬ 
marized as follows g] . Suppose an atomic system with Ni 
ground levels and N 2 excited levels (N = Ni + N 2 ) reso¬ 
nantly interacts with a laser held. Under an assumption 
of no damping, the interaction Hamiltonian is written in 
the rotating-wave approximation (RWA) gHEl as 

H - (y, o) • « 


where V is an N\ x N 2 coupling matrix between the 
ground and excited levels. Then, the MS transformation 
theory predicts that the coupling matrix V is diagonal¬ 


ized through a MS transformation given by 

H' = UHU (2) 
where the unitary transformation U defined by 

u = (* “) . (3) 

Here the block matrices X and Y are the unitary trans¬ 
formations that diagonalize VV 1r and V^V, respectively. 
The resulting Hamiltonian H' is a direct sum of M sets 
of two-state systems and N — 2 M uncoupled single states 
(2 M<{N u N 2 }), i.e., 


H' 




f o n 2 \ 
[n * 2 o J 


( 0 

Wm 0 J 


® (hi) ® (h 2 ) ® • • ■ ® (hpj— 2 M ), (4) 


where $Vs are the Rabi frequencies, that are all different 
in general, of the two-state systems. 

Various systems have been analyzed by the MS trans¬ 
formation [2 ] , which include two-state superposition sys¬ 
tems three-state A-linkage systems 0, four-state 

diamond- and tripod-linkage systems [TOlfM] . V-state 
M- and TV-linkage systems g], to list a few. The 
MS transformation of three-state A-linkage systems is 
in particular interesting in regards to the laser tech¬ 
nique called stimulated Raman adiabatic passage (STI- 
RAP) [ig, where the systems undergo an adiabatic evo¬ 
lution between MS-transformed coupled two states, leav¬ 
ing the dark state never populated during the laser in¬ 
teraction. The result is a complete population transfer 
(cpt) gmnni] from one ground state to the other, 
which is easy to understand in the context of the MS 
transformation. 

In this paper, we describe an experimental study of 
MS-transformed IV-state systems performed by ultrafast 
laser-induced coherent excitations. The systems under 
consideration are the complex coupling linkages of the 
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D1 and D2 transitions of atomic rubidium ( 85 Rb). For 
example, the D1 transitions from |5 Si/ 2 ,F = 2,3) to 
\§Pi/ 2 ,F' = 2,3) correspond to an N = 24 system with 
N\ = 12 and N 2 = 12. We first analyze the atomic 
rubidium systems in terms of the MS transformation to 
derive the Rabi frequencies for the D1 and D2 transitions 
induced by an elliptically-polarized light in Sec. II. Af¬ 
ter briefly describing the experimental procedure in Sec. 
Ill, it is shown based on our experimental results that 
the system dynamics are all governed by a sum of at 
most two decoupled Rabi oscillations [18) [19] when these 
systems interact with ultrafast laser pulses of any polar¬ 
ization state, in Sec. IV. The implication of this result 
to possible controls of photo-electron polarizations is dis¬ 
cussed before the conclusion in Sec. V. 


II. THEORETICAL CONSIDERATION 

A. Linear and circular polarizations 

We first consider a laser pulse in the linear or circular 
polarization state interacts with the atomic rubidium. 
Suppose the polarized electric field is given by 

E(t) = e£(f) cos(wf + </>), (5) 

where the polarization vector e = z or (x±iy)/y /2 for the 
linear or circular polarization, respectively. The D1 tran¬ 
sitions is dictated by the selection rule of Am = 0 or ±1 
(q = Am ), so, for each m, the Hamiltonian can be writ¬ 
ten in the four-state basis {|F = 2,m),\F = 3,m), |F' = 
2 , to ± g), | F' = 3, to ± q}} as 

/ 0 0 fJ-22'E(t) p.23 'E(t)\ 

rr( + \ _ 0 & P32'E(t) M33 ’E{t) I 

M22 >m* M32 W 0 O ’ 

M33 >E(t)* 0 WO + S'J 

( 6 ) 


where 8 — luf =3 — top— 2 , S' — ujp/ = 3 — ujpr— 2 , cuq — 
ojf ’ = 2 — u>f = 2 , and the dipole moments are given by 
Hff’ = (F,m\er q \F',m-q). 

Since we consider ultrafast laser interactions, the in¬ 
teraction time is extremely short compared to the in¬ 
verse of any hyperfine energy splitting, so the effective 
Hamiltonian under a resonant approximation [20] can be 
alternatively written in the atomic frame as 


H{t) 


/ 0 0 A B\ 

0 0 CD 

A* C* 0 O ’ 
\B* D* 0 0/ 


(7) 


where A = /i 22 '£(t)/2, B = Ai 23 '£(t)/2, C = /r 32 /£(t)/2, 
and D = /X 33 ' £(t)/2. When <f> = 0 is assumed without 
loss of generality, the coupling matrix V ( t ) is given by 


v W ~ (c d) ’ ^ 

and the MS transformation decouples the system into two 
independent two-state systems through diagonalizing the 
matrices V(t)V^(t) and V^(t)V(t) that are given by 


v{t)vH0 = v\t)v(t) = + f D + b d ?). 

(!>i 

Then, the coupling strength of each two-state system, or 
the square of each Rabi frequency, is determined by the 
eigenvalues | 2 ] given by 


A 2 + B 2 +C 2 +D 2 ± v /( A 2 + B 2 + C 2 + D 2 ) 2 - 4 (AD - BC) 2 

^- ■ ( 10 ) 


In the basis of the eigenvectors |A±), an arbitrary ground 
state can be a pure superposition state or a mixed state 
of |A±). If these two eigenvalues A± are the same, the 
two independent oscillators represented by |A±) vibrate 
coherently with the same Rabi frequency, thus the CPT 
could occur. Here, the condition for A+ = A_ is given by 

{{A + D) 2 + {B - C) 2 }{{A- D) 2 + (B + C) 2 } = 0. (11) 

As for the n transitions {q = 0), Clebsch-Gordon coef¬ 
ficient symmetries [23] allow ( 2 TO| 2 | 2 m') = —(3m|z|3?n') 


(A = —D) and (2?n|z|3m') = (3m\z\2m') (B = C ) for all 
m = m' , which satisfies A+ = A_ and the Rabi frequen¬ 
cies are the same for all m’s. Therefore, the CPT of the 
total system occurs and, moreover, the system initially in 
any coherent-superposition or mixed ground-state sim¬ 
ply undergoes simply a Rabi oscillation with the Rabi 
frequency given by 

n(t) = yjn\ 2 , + £.(€)/2K. (12) 

The Rabi oscillations in these complex degenerate-level 
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manifolds have been observed in an implementation of 
single qubit operations of atomic rubidium |2lj . 

Meanwhile, for the a transitions (q = ±1), the equal¬ 
ity AD = BC holds, which makes one of the eigenvalues 
zero and the other twice the eigenvalue of q = 0 (the tt 
transition). In other words, the a transitions for differ¬ 
ent to’s have respectively a pair of uncoupled states and 
their Rabi frequencies are y/2 times bigger than that of 
the corresponding 7r transition. Therefore, the temporal 
dynamics of transitions depend on the initial con¬ 
dition, unlike the 7r transitions, since its oscillators are 
drastically different, which behavior has been applied to 
coherent controls of medium gains [22] ■ 


B. Elliptic polarizations 

Now we generalize the problem by considering ellipti- 
cally polarizations of the laser interaction |3). When the 
electric held is defined by 

E(f) = d+£+(t) cos (ut + <j)+) + cr_£_(f) cos (ut + (/_), 

( 13 ) 

the elements in V are given by £±(t)r]™ e l< ^ ± /2, where 
are the transition dipole moment, ( Fm\er±i\F'm'). 
The dipole moments can be factorized by a common fac¬ 
tor, [ijji = (J = l/2||er||J' = 1/2) for all {F,m} (23] , 
When the ellipticity e is defined by (£+ —£?_)/(£+ + £;!), 
the Rabi frequencies are given by 


n± = if2 0 


lie 


(14) 


with Q 0 = /rjj'(£^_ + t 2 _Y^ 2 / 2 h. 

When the Hamiltonian is MS-transformed through 
the summarized procedure in Eqs. (l)-(4), the resulting 
Hamiltonian H' is a direct sum of twelve sets of two-state 
systems given by 


H’ 


( 0 n +e ^+\® 6 

\n+e-**+ 0 J 

f 0 

0 ) ’ 


(15) 


III. EXPERIMENTAL PROCEDURE 


The experimental setup is schematically shown in 
Fig. [l] Rubidium atoms ( 85 Rb) were cooled and trapped 
in a conventional magneto-optical trap [28]. The detail 
of the experimental setup was reported elsewhere m 
21 [33. ■ Then the atoms were interacted with ultra- 
short laser pulses from a Tksapphire mode-locked laser 
amplifier operated at a pulse repetition rate of 1 kHz. 
The center wavelength of the laser pulses (the pump 
pulse) was tuned to 794.7 nm (or 780 nm) for the rubid¬ 
ium D1 (D2) transition with a full-width-half-maximum 
(FHWM) bandwidth of 3 nm (0.5 nm). The single pulse 
energy was varied up to 0.03 mJ, which was equivalent 
to a Rabi pulse area of 0 = 27 t when the pulse diameter 
at the atom was about 0.5 mm and the size of the atom 
cloud was 200 /im. The polarization was controlled by a 
broadband quarter-wave plate (QWP). 
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FIG. 1: (Color online) (a) Schematic diagram of the exper¬ 
imental setup, (b) Magneto-optical trapping apparatus for 
85 Rb. (c) Energy level diagram of 85 Rb atoms. 


where the first set of six two-state systems corresponds 
to the a + transitions and the second to the er~ transi¬ 
tions. This result can be readily understood in terms of 
the basis reduction from |F, tof) to J, rrij) [24]. Since 
spin interactions (hyperfine interactions in this case) are 
irrelevant in ultrafast-time scale dynamics, they can be 
traced out, and the reduced basis description in J, nij) 
can be alternatively acquired. The advantage for using 
the MS-transformed \F,mp) bases, instead of using the 
reduced basis set, comes from the fact that the coherent 
dynamics coupled with spin interactions can be directly 
seen in the former bases, implying applications to ultra¬ 
fast time-scale nuclear-spin polarization controls mm- 


The measurement of the excited state (5 P 1 / 2 or 5 P 3 / 2 ) 
population of the atoms was achieved by photo-ionizing 
the excited atoms by frequency-doubled ultrafast pulses 
(the probe pulse) from the same laser. The pump and 
probe pulses were independently focused and delivered 
by a dichroic mirror (DM) to the MOT. The number of 
ions (Rb + ) was counted by a micro-channel plate detec¬ 
tor (MCP). The signal linearity was assured by operating 
the experiment in the one-photon perturbation regime. 
The error due to three-photon ionizations by the pump 
pulse was estimated below 5% at a pulse-area of 37 t (30], 
which was negligible in the considered experimental con- 
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dition. The whole experiment (four steps: MOT-turn off, 
laser-control, ionization, MOT-turn on) was repeated at 
2 Hz by cyclically turning on and off three mechanical 
shutters. 


IV. RESULTS AND DISCUSSION 

The MS transformation of the rubidium D1 transitions, 
described in Sec. |TI] predicts the relation between the 
Rabi oscillations of the 7r and a transitions given, in terms 
of the Rabi-oscillation pulse-areas, by 

©7T = (16) 

which relation is comes from Eq. ( fl4| that gives f2(e = 
0) = n o /V3 and H(e = ±1) = ±^/2/3n o or zero. Also, 
with the same reason, a half of the ground levels remain 
intact for the a transitions if the atom is unpolarized, 
resulting in the relation of the maximal excitation prob¬ 
abilities given by 


= 2P“ 


(17) 


To verify these predictions, Rabi oscillation experiments 
of rubidium D1 transitions performed by 7r- and a- 
polarized ultrafast laser pulses, respectively, were per¬ 
formed. The result is shown in Fig. [2j clearly exhibiting 
that the Rabi frequency of the a transition is \/2 times 
larger than that of the the 7r transition and that the os¬ 
cillation amplitude of the a transition is a half of that of 
the the 7r transition, thus the two relations in Eqs. (16) 
and are confirmed. The theoretical calculations are 
respectively plotted with dashed lines, of which the devi¬ 
ation from the experimental result is due to the spatially 
inhomogeneous laser interaction with the atom ensem¬ 
ble 30]. When the spatial averaging effect is taken into 
account (by assuming the Gaussian laser beam and atom 
ensemble size ratio is 2.5), the new calculations, shown 
with solid lines, agree well with the experimental results. 


In the second experiment, elliptical polarizations are 
considered. When the rubidium atoms interact with an 
ultrafast laser pulse in an elliptical polarization, the inter¬ 
action Hamiltonian in Eq. (J7| is expressed as a superpo¬ 
sition of cr + and <j~ transitions. In this case, the decou¬ 
pling is not simply attained through a 4-by-4 MS trans¬ 
formation H' = UHW in Eq. (|3]) . However, we can still 
describe the dynamics of the MS-transformed two-state 
system as a mixture of the Rabi oscillations that corre¬ 
spond to the constituent polarization-specific transitions. 
In the D1 transitions performed by an elliptically polar¬ 
ized, there are parts of the two-state system that are re¬ 
spectively run by a + and cr” transitions. In that regards, 


we can decompose the initial ground-state atom, writ¬ 
ten in | J,mj) bases, into | ip±) = a± |l/2, ±1/2), where 
|a + | 2 ± |a_| 2 = 1, | i/j±) are the Rabi oscillation bases, 



Pulse area @ CT 


FIG. 2: (Color online) Ultrafast Rabi oscillations of 85 Rb 
atoms between 55 'i/ 2 and 5 Pi/ 2 energy states. Experimen¬ 
tal results for linearly (the 7r transition) and circularly (the 
a transition) polarized ultrafast laser pulses are shown with 
blue open circles and red filled circles, respectively. The two 
results are commonly scaled along the vertical axis by the 
maximal population of the 7r transition and along the hori¬ 
zontal axis by the 7r-transition pulse area. Black dotted lines 
are the theoretical calculation based on Eq. (181 for the cor¬ 
responding Rabi oscillations. New calculations that take the 
spatial average effect into account are represented in solid 
lines. 


respectively, for transitions. Therefore, the dynam¬ 
ics of the excited-state population of the atom under an 
elliptically polarized light is given by 


Pdi = |a+| 2 sin 2 © + + |a_| 2 sin 2 0_, 


(18) 


where 0± = a/( 1 ± e)/30 are the parametrized 

pulse areas responsible for the o p± transitions, re¬ 
spectively, obtained from Eq. © , with 0 = 

(2fi,) -1 (1/2 ||er|| 1/2) £(t)dt. The tt transition is 

naturally defined for e = 0. Figure [3] represents the 
experimental results. The excited state population in 
5 Pi/ 2 appears as a degraded oscillation in each ellipticity 
e £ {0.3,0.4,0.45} as shown in Figs. [3](a-c) , because, as 
in Eq. (18), it is given by a sum of two distinct oscil¬ 
lations. The calculation based on Eq. and consid¬ 
ering inhomogeneous pulse areas for various ellipticities 
is represented in Fig. [3](d) . The Fig. [3])e) represents the 
experiment result that close up the first peak for various 
ellipticity. When the ellipticity is varied from zero to one, 
the first oscillation appears at a gradually smaller pulse 
area, and the pulse area for the first peak is changed from 
7T to tt/v^, as shown in Fig. Did, e). 
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FIG. 3: (Color online) Experiment results for D1 transitions by elliptically polarized lights, (a-c) The blue circled lines represent 
experimental data for (a) t = 0.3, (b)e = 0.4, and (c)e = 0.45 . The black dotted lines are the ideal Rabi oscillation curves 
and the blue solid lines are their spatially-averaged oscillations, (d) The calculated Rabi oscillations for various ellipticity are 
represented in a three-dimensional plot, (e) The experimental data of first peak position are closed up for various ellipticities. 


Lastly, Rabi oscillations of D2 transitions between 
5 Si/ 2 and 5 P 3 / 2 energy states are considered. In 
Fig. [4]) a), the experimental results are compared with 
theoretical calculations. After a similar argument lead¬ 
ing to Eq. ©, the excited state population as a result 
of D2 transitions is also given by a sum of two Rabi os¬ 
cillations as 

P D 2 = l a +| 2 sin 2 \j~Y~® + |a-1 2 sin 2 \j =-^0, (19) 

where 0 = (2fi) -1 J'^ oc (1/2 ||er|| 3/2) £(t)dt and the 
a± denote the initial probability amplitudes of the 
ground state in | J = 1/2, mj = ±1/2) bases. The 
Clebsch-Gordan coefficients for their coupling to 
| J' = 3/2, m'j = ±) bases are shown in Fig. B b )- 
However, a direct comparison of the excited-atom popu¬ 
lations in 5P 3 /2 levels by counting photo-ions results in 
an experimental artifact, because the ionization passages 
via one-photon transitions to the continuum, as shown in 
Fig.[4|c), for | J' = 3/2, m'j = ±1/2) and | J' = 3/2, m'j = 
±3/2), are different from each other. The ratio be¬ 


tween the ionization rates may be estimated by compar¬ 
ing the bound-to-bound transitions of the corresponding 
angular symmetries. Since the coupling strengths from 
| J = 3/2, mj = 1/2) and | J = 3/2, mj = 3/2) are pro¬ 
portional to J2j> 1(3/2 \\ez\\ kJ ') ( kJ ', l/2|3/2,1,1/2, 0)| 2 
and Yhj' 1(3/2 ||er||/cJ , )(fcJ',3/2|3/2,1,3/2, 0)| 2 , respec¬ 
tively, when |(J||er|| kJ ')\ 2 oc (2 J' + 1)/(2J + 1) is as¬ 
sumed, it can be estimated that the ratio of the ionization 
rates is 1.06, which is in a good agreement with 1.1 used 
to fit the data to the theoretical calculations in Fig. |4])a). 

The remaining experimental errors are due to laser 
power fluctuations and the center mismatch between 
laser and atoms, both of which are predominant devi¬ 
ations at large pulse areas. The laser fluctuation within 
a 10% of a shot-to-shot deviation was kept low by event 
statistics. The size of the atom cloud (200 /jm) was small 
compared to the steering mirrors distance (300 mm), so 
the intrinsic misalignment of 30 /im due to steering an¬ 
gular resolution led to an experimental error. However, 
the guidelines that considered the spatial averaging effect 
are valid in the experiment m ■ 
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FIG. 4: (Color online) (a) Experiment results for D2 transi¬ 
tions. Experimental data are represented by blue open circles 
for circular polarizations and red closed circles for the lin¬ 
ear polarization. In comparison, theoretical calculations and 
their spatially-averaged results are shown with dotted lines 
and solid lines, respectively, (b) Coupling strengths defined 
in | J, mj) bases, (c) Schematic comparison between the ion¬ 
ization schemes for 5Pi/2 and 5 P 3/2 states. 

We now turn our attention to the implication of the 
obtained results to possible applications towards polar¬ 
ization controls of electron or nuclear spins. Electron 
spin polarization controls by photo-excitation of atoms 
have been theoretically discussed for alkali atoms (31] . 


where coherent excitations to nP states by circularly- 
polarized light was considered so that spin-orbit coupling 
exchanged spin and angular momenta. Alternatively, it 
can be considered, based on our demonstration, that the 
excited nP state is independently controlled to nPi/ 2 and 
nP 3 / 2 through ultrafast Rabi oscillations. So, for exam¬ 
ple, the procedure to control the electron polarization 
can be as follows: First, a circularly-polarized pulse res¬ 
onantly tuned to nPi / 2 state completely excites the pop¬ 
ulation in a particular electron spin state, for example 
|nij = 1/2), leaving the rest intact. The second pulse 
resonantly tuned to nP 3 / 2 state completely excites the 
remaining ground-state atoms, which are then photo- 
ionized by an extra means of an ultrafast time resolu¬ 
tion. Finally the third pulse identical to the first pulse 
de-excites the atoms in nPi/ 2 back to the ground state 
to complete the spin polarization control. Then, it is 
achieved that every ionized photo-electron has pure spin 
states. Higher fidelity controls may be achieved through 
adiabatic passages by applying frequency chirps. 

V. CONCLUSIONS 

In summary, we have studied ultrafast-laser driven 
Rabi oscillations of MS transformed atomic rubidium 
transitions. The complex coupling linkages involved with 
the D1 or D2 transitions are analyzed based on the 
MS transformation and the simplified decoupled systems 
are experimentally investigated by ultrafast laser interac¬ 
tions of spatially-localized rubidium atoms in a MOT in 
the strong interaction regime. It has been found that all 
the system dynamics are governed by a sum of at most 
two decoupled Rabi oscillations for any polarization state 
of the interacting ultrafast pulse. The implication of this 
result has been discussed to control photo-electron polar¬ 
izations. 
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